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ABSTRACT 

Two seemingly unrelated scientific disciplines, information processing and 
quantum mechanics, were separately developed, until physicist Richard 
Feynman proposed their combination. Today, quantum computation is at the 
forefront of research in theoretical physics and computer science. 

The objectives of this thesis are to present an introduction to the subject 
of quantum computation and its potential, and to analyze the physics underlying 
the operation of one of the proposed realizations of a quantum logic gate, based 
on cavity quantum electrodynamics (QED). We discuss the fundamental 
aspects of quantum mechanics that may be a revolutionary new way to process 
information. We then present the appropriate theoretical background of QED 
needed to understand the proposed realization of a quantum gate. We identify 
the advantages and disadvantages of this realization. We then address and 
discuss the phenomenon of decoherence, the most fundamental obstacle to 
overcome for the development of practical quantum computers. Finally, some of 
the possible military applications are presented, along with our thoughts for the 
future of this field. 

We conclude that quantum gates based on cavity QED are feasible with 
current technology. However, the construction of a quantum network with these 


gates is not feasible, mainly due to the decoherence problems. 
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I. INTRODUCTION 

Every two years for the past fifty, computers have become twice as fast 
while their components have become twice as small. The main reasons for this 
progress are the technological advances made in lithographic manufacturing 
techniques. The implementation of these techniques allows us to build sub- 
micron wide logic gates and wires, cramming them onto the surface of silicon 
chips. However, integrated circuit technology is running up against fundamental 
limits, due to the laws of physics. Even if we find new ways to manufacture 
smaller chips, ultimately they will be so small that they will be made out of a few 
atoms. Reaching the atomic scale, chip manufacturers will face the rules of 
quantum mechanics, which are quite different from the classical rules that 
determine the properties of today’s conventional logic gates. So if computers are 
to become smailer in the future, it is inevitable that quantum mechanics will play 
an important role. 

The point is, however, that quantum theory can offer not only the means 
for new technological breakthroughs in construction techniques, but it can 
support an entirely new kind of computation, the so-called quantum computation. 
On the other hand, the combination of quantum mechanics and computer theory 
can offer a new way of understanding quantum systems themselves. 

That was the starting point for this new theory, as the physicist Richard 


Feynman proposed in 1982 a scheme for the simulation of quantum-mechanical 











systems by other quantum systems, built for this purpose.” The next step that 
revealed the power of quantum computation was in 1985 when David Deutsch of 
the University of Oxford published a theoretical paper® in which he described a 
universal quantum computer and its potential. Simultaneously Deutsch and other 
physicists realized that quantum mechanics offered a way of overcoming one of 
the most fundamental difficulties of computer science, that of a unique 
computational complexity for every mathematical problem.* It was believed until 
then that the time needed to solve any given problem that was polynomial or 
greater in the size of its inputs was independent of the physical apparatus used 
to perform the computations. This axiom seems not to be true for quantum 
computers. It was this realization that made scientists all over the world start 
looking at the development of these new ideas. 

For more than a decade after the first proposal the obstacles encountered 
were so many and difficult that it seemed this new theory would be yet another 
academic curiosity, without ever offering any of the initial believed benefits. 

But in 1994, Peter Shor from AT&T's Bell Laboratories in New Jersey 
devised the first quantum algorithm that, in principle, could perform efficient 
factorization’, something useful that only a quantum computer could do. Today, 
most of the methods used for encryption are based on the difficulty of 
factorization, thus allowing the privacy and security in our electronic transactions 
(bank accounts, credit cards, etc.). The fact that a quantum algorithm could 


solve this problem gave a new boost for the prospects of quantum computation. 





Since then a growing number of scientists and institutes all over the world 
have developed various projects trying to explore the theoretical background 
and the technology needed to build a quantum computer. 

In this thesis, first we analyze the elements of quantum mechanics that 
give a new way to process information and compare it wherever applicable with 
classical computational methods. Then we will briefly discuss some of the 
proposed realizations and applications for quantum computers. From these 
proposals we will focus on the one that uses the principles of cavity quantum 
electrodynamics to build the necessary quantum logic gates. We will review in 
Chapter IV the theoretical background needed for the implementation of these 
techniques. Then in Chapter V we will report some of the recent experiments 
based on cavity quantum electrodynamic techniques and we will describe a 
proposed 2-bit quantum universal gate. In Chapter VI we will discuss the most 
difficult obstacle needed to overcome for building a quantum computer, the 
phenomenon of decoherence. Finally in Chapter VIl we present some possible 
military applications, while in Chapter VIII we discuss the possibility of building a 
quantum computer in the near future. 

It should be noted that quantum computation is a part of a new evolving 
science called quantum information, which also includes quantum teleportation 
and quantum cryptography. Even though some of the aspects of the other two 
fields can be combined with quantum computation or can be used to solve some 


of the encountered problems, it is beyond the scope of this thesis to provide a 





review or analysis for these fields. In some cases there may be a review of the 


advancements in these theories, which provide a deeper understanding of some 


of the aspects of quantum computation. 





Hl. ELEMENTS OF QUANTUM COMPUTERS 


Quantum computers are different from conventional ones due to their 
distinguished fundamental elements and their rather peculiar properties, which 
arise from quantum theory. 

The most important elements of quantum mechanics that give these 
properties to quantum computers are the superposition of states, quantum 
entanglement, reversible unitary evolution, and decoherence. We will review 
these and explain their effects in the construction of a quantum machine. A 
quantum computer is made by quantum logic gates, which operate with quantum 


bits. The gates are connected with quantum “wires”. 


A. QUANTUM BITS AND REGISTERS 
1. Definitions and Properties of Qubits and Registers 


The quantum bit or simply called “qubit” is.the building block of quantum 
logic, like the bit, which is the building block of conventional logic. Any quantum 
system that can have two accessible states can be used, at least theoretically, 
as a qubit. However, the property of superposition distinguishes a qubit from a 
classical bit. A classical bit can have only two choices 0 or 1 and is always in 
either of these states. The qubit, however can have an infinite number of 
superimposed states between the assigned values of 0 and 1. Its properties are 
based on these infinite states. Using mathematical terms we can say that the 


state of a quantum system is a vector in Hilbert space that can be written as a 











superposition of basis vectors, chosen by us with certain mathematical 
properties.° For a two-state system in which we choose as basis vectors the 
classical realization for the states 0 and 1 we can write: 

|v) =a|0)+|1) (1) 
where a and b are normalized complex coefficients. Due to this definition the 


value of the qubit is uncertain. The probability to find the qubit in state O or in 
state 1 is given by |a|’ or |4|’ respectively. It should be noted here, that 


according to quantum theory, the measurement problem of a quantum system 
imposes a limitation on the accessible amount of information. Most of the 
information included in the superimposed states is not accessible to us, but it is 
nevertheless needed in order to predict the evolution of the system until its final 
state, which includes the solution.’ 

Since the qubit has a basis of two states we can infer that n qubits will 
have a basis of 2" states. The collection of n qubits is called a quantum register 
in direct analog with a classical computer register, which is a collection of n bits. 
A classical register has one definite value even though it involves n_ bits. 
However, a quantum register consisting of n qubits can exist in a superposition 
of all possible classical states, thus allowing the performance of a single 
computation to solve a problem with many different inputs. This property is 
called quantum parallelism and gives the advantage of solving a problem in less 


time by utilizing a quantum machine.® 





2 Entanglement 


The notion of superposition for the n qubit register described before was 
very simple and did not include a fundamental property of quantum theory, 
called quantum entanglement. 

Since 1935 wien Einstein, Podolsky, and Rosen presented the famous 
EPR paradox’, an endless debate has started considering entangled particles 
and their properties. 

Schrédinger first used the term entanglement in one of his papers” in 
1935, where he commented on the EPR paradox.® Two equivalent definitions for 
entanglement®"' are: 

The states of the subsystems constituting a complete quantum system do 
not determine the state of the system. 

Some definite states of a complete quantum system do not correspond to 
definite states of the subsystems. 

We will try to describe entanglement for a system consisting of two 
subsystems®*. Each subsystem has its own Hilbert space denoted by Ha and Hg. 


We assume that the vectors |), and |m)_ are a complete orthonormal basis for 


the Hilbert spaces Ha and Hg respectively (with n,m=0,1,2,..). The total Hilbert 
space for the complete system is given by the tensor product of the individual 


Hilbert spaces of each subsystem, 7, ®H, or |n), ®|m), =|n) |m),. So now we 


can write any state of the complete system as: 











[Y) se = DSnm|?2) 417), (2) 


where Cam are normalized complex coefficients with Yleam| =1. Now we define 


nm 


the direct product state as a state |v) ie of the complete system that can be 


factored as a tensor product of two normalized states |y), = >e,|n), and 


lv), = Lene 


¥) 0 =1¥)IY)e (3) 

We can see that there are states in the Hilbert space H,®H, that are 
not direct product states. These are called entangled states. An example is the 
state |0),|0), +[1),[})s- 

The property of entanglement is encountered in every aspect of quantum 


information theory and has been used in many applications such as quantum 


coding, quantum compression of information and quantum teleportation. 


B. QUANTUM GATES 


Logic gates perform information processing in a classical computer. The 
logic gates change the input bit in accordance with a truth table, defined by 
binary logic. 

The quantum gate is implemented by means of a unitary evolution 


operator that changes the input state in accordance with the same truth table as 


in classical gates.°" 








Classical gates perform operations on definite states and have definite 
outputs, while quantum gates can perform operations in a superposition of input 
states resulting in a superposition of output states. This is the notion of quantum 
parallelism described before. 

The key feature that makes quantum gates powerful is the unitary 
evolution, which ensures the isolated performance of operations and their 
reversibility. However, this feature is the most difficult obstacle to overcome in 
the implementation of quantum computing, since it is not easy to prepare a 
quantum system and make it evolve without any interference from the 
environment. This difficulty will be discussed in more detail later. Now we will 
analyze the notion of this unitary evolution operator and its correspondence with 
a quantum gate. 

1. Reversible Unitary Evolution Operator 

By definition® an operator U is unitary if its inverse U' is equal to its 
adjoint U*, i.e. U’U=UU"=I. If H is a Hermitian operator then the operator U=e" is 
unitary. 

An isolated quantum system evolves in a reversible manner, according to 


the time dependent Schrédinger equation: 
ga 
ih |v) = HOlv) (4) 


Since equation (4) is linear and homogeneous we can infer that the 


correspondence between two states at different times is linear. Therefore we can 








introduce a linear operator U(t, to), which evolves a state at time tp to a state ata 
later time t: 
ly) =U(6t.)|v@)) (5) 


If we substitute (5) into (4) we have: 
ihoU( 1, )w(to)) =H QUE 4)|¥ (to) => ihoU( ,)=HOUGt) (6) 


When the operator H does not depend on time we can obtain by 
integration of equation (6): 

U(t,t.) =e" (7) 
and so the state at a time t can be written: 
lyn) =e" w(t,)) (8) 

The operator U is unitary and a quantum gate will perform calculations 
according to this operator, thus transforming the input state at time to to an 
output state at time t. The operator H is the Hamiltonian of the closed quantum 
system. It does not involve any effects from the environment. It is now clear the 
reason why the system must be closed. Also it must be noted that while the 
computation is performed we can never know the states the gate runs through 
until we make the measurement at a time t. This has the consequence that a 
computation can be performed by many different ways, even though we always 


achieve the same final result’. 


10 








Another fine point involving the unitary evolution operator is the 
dependence or not of the Hamiltonian H on time. Above we assumed that H was 
independent of time and we derived equation (8) for the evolution of the system. 
However, this is not the case when H depends on time. A lot of books discussing 


quantum computing define the evolution operator in this case as: 


_t 
- f H(t')dt’ 


U@t,)=e ° 


FQ) 


which is not always true®. The derivative of an operator of the form e'™ is equal 


to F’(t)e*only when F(t) and F(t) commute. So if someone chooses to 
implement a quantum gate with a system described by a time dependent 
Hamiltonian, they must ensure the above commutation. When the Hamiltonian 
does not commute with itself at different times, the time evolution operator may 
be defined in terms of the Dyson time-ordering operator. ° 

The evolution of a closed quantum system described by a unitary 
operator provides us the theoretical background for the operation of a quantum 
logic gate. A system with two accessible states can be used, at least 
theoretically, to implement such a gate.* Many two state systems are known in 
physics. An example is the spin up and spin down states of a spin-1/2 
elementary particle. 


2. Density Operator 


It often occurs that the state of a system is not perfectly determined or 


where a large number of identical quantum systems are used for the 


11 








implementation of a quantum gate. In these cases the evolution of the system 
can not be described by the unitary operator defined in the previous paragraph. 
However, there is a mathematical tool called the density operator, which ensures 
the simultaneous application of quantum mechanics and probability theory. The 
density operator is an alternative way to describe any quantum system. It also 
provides the basis for understanding the irreversibility problems involved in the 
evolution of quantum systems.” We first introduce the notion of density operator 
for a system whose state is perfectly known. This is called the pure state.° Then 
we will expand the definition for the non-pure case, where we have a statistical 
mixture of states.° It must be noted that in our treatment of various elements of 


quantum mechanics we will use the Dirac notation. 


a) Density Operator in the Pure Case’ 


We consider a system with state vector: 


lv@)= de. Ole.) (9) 





where |w,) are a complete set of orthonormal vectors in Hilbert space and the 
coefficients satisfy the relationship: 


Yle,@) =1 (10) 


The matrix elements of an observable A can be written as: 


and its mean value is given by: 


12 








(A) = (¥@|AlvO) = 216, On 4m (12) 


Equation (12) shows us that mean values of observables are given 


as quadratic expressions of c;(f)c,,(f). But from the definition of the state vector 


we can write: 





CF (De, (t) = (Un WOXYO|4,) (13) 


Now we define the density operator p(f) as: 
ep) =|yOKwO| (14) 
This can also be represented by a matrix called the density matrix, 


with elements: 





Prn(t) = (Um |PO)|4n) = Or On) (15) 
Now that we have defined the operator we can find how a quantum 


state of the system is characterized. Using equation (15) equation (10) becomes: 
DYle.O! = DT Pm = Tre@ =1 (16) 


where 7r[p(1)]is the trace of the matrix. 
Using equations (11) and (15) equation (12) becomes: 


(AD) = >i (u, PO u,\u, \Alu,,) = >. (4, POA#,,) = HPO] (17) 


nm m 








Finally if we use Schrédinger’s equation we can write the evolution 


of the operator as: 


in pO =HO.PO| (18) 
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The density operator is hermitian and for a pure state it satisfies: 


p= pl), Tre’ O=1. (19) 


The last two relations are true only for the pure case. 


b) Density Operator for a Statistical Mixture of States (Non- 
Pure Case) 


Now we consider a system whose state is not perfectly known. In 


this case we assign to each probable state of the system |v.) a probability px, 
with >) p, =1. 
k 
We define the density operator for the state lv.) as p, =|Vi Wel. 
Then the density operator at any time t is: 


pit)= > PP) (20) 


The trace of this operator is: 


Trp =>, Trp, = DP: =) (21) 
k ke 


The evolution of the operator is again given by equation (18) 
above. It is noted that the last two relations in (19) do not hold and we have: 
p’#p, rp’ <1 (22) 
The density operator and its properties will be used in Chapter VI 


when the phenomenon of decoherence will be presented. 
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3. Universal Quantum Gates 


As noted earlier one of the most important features of aentint gates is 
their reversibility. David Deutsch has defined’* the universal quantum logic gate 
as the gate that can be used to simulate any quantum gate. Since early 1980s a 
lot of proposals have been published for such gates. It has also been shown that 
for classical reversible computation the simplest universal gate has three inputs 
and three outputs. From all the proposed schemes two gates have received the 
most attention as potential candidates for implementation in quantum computing, 
the Toffoli® and Fredkin’® gate. 

In the Toffoli gate two input bits control the state of the third bit (target). 
The control bits do not change between input and output. The state of the target 
bit changes when the two control bits are set to 1. 

In the Fredkin gate there is one control bit and two target bits. The control 
bit does not change. The target bits swap states only when the control bit is set 
to 1. Tables 1 and 2 are the truth tables for the Toffoli and Fredkin gates, 


respectively. 





Table 2.1 Toffoli gate 
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Table 2.2 Fredkin gate 


Even though a variety of systems have been proposed for the 
implementation of quantum computation only a few of them can exhibit their 
usage as a universal quantum gate. The difficulty of implementing systems with 
three bits has been reduced since 1995, when a team of scientists (Tycho 
Sleator and Harald Weinfurter) described a 2-bit universal quantum gate.” A 
2-bit gate can be used as the building block for any quantum logic network. This 
was shown at the same time by other researchers as well (i.e. David di 
Vincenzo). The proposed scheme for the implementation of this gate is based on 
cavity quantum electrodynamics theory and specifically to the interaction of 
atoms with single mode fields of microcavities. 

We have not yet mentioned a way for connecting the quantum gates. The 
so-called quantum wires are an important part for the construction of a quantum 


computer. Some proposals exist for such wires*, but we are not going to 


describe any of them in this thesis. 
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Hl. PROPOSED APPLICATIONS AND REALIZATIONS FOR QUANTUM 
COMPUTERS 


We present in this chapter some of the proposed applications for 
quantum computers. We also present a brief description of the basic proposed 
schemes for the implementation of quantum computers. Since the power of 
quantum computing is now starting to be revealed there will surely be more 
proposed applications in the future. Some of our thoughts for possible military 


applications will be presented later. 


A. APPLICATIONS 


Current conventional machines are capable of performing very complex 
computations and solving many difficult problems. However, there are some 
problems that either appear to be intractable or the time and computing power 
needed for the solution is enormous. Thus, most of the researchers in the field 
focused their efforts in quantum computers that could solve these special 
problems. 


1. Factorization 


As was noted earlier the difficulty of factoring large numbers is the key 
issue for efficient cryptography. Even though the factorization of a large integer 
cannot be considered classically intractable, it is very difficult since it requires 
huge computing power and a very large number of calculations. This is due to 


the insufficient classical algorithms for the solution of the problem. 
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In 1994 Peter Shor from AT&T's Bell Laboratories in New Jersey 
devised the first quantum algorithm that, in principle, could perform efficient 
factorization of large numbers.° This event was a major breakthrough and the 


factorization problem became one of the most interesting applications for 


quantum computers. 


Shor’s algorithm is based on the fact that factoring a number can be 
related to the problem of evaluating the period of a function.® Then the factors of 
the number can be found using number theory. The evaluation of a period of a 
function by classical computers is very difficult, since it requires the same 
computational power as the factorization problem. Quantum computers can 
solve this problem rather easily by using their advantage to perform the required 
calculations in parallel. The application of a discrete Fourier transform in the 
final state of the quantum system yields the period of the function. 


2. Database Search 


Another problem that can be solved by quantum algorithms faster than 
conventional ones is the search of a database. If someone wants to find a 
particular item searching a random list of N items has to examine at least N/2 of 


them in order to ensure a 50% success probability.° In 1996 Lov Grover 


showed" that a quantum computer could perform this task in VN steps. This is 


an increase of speed by a factor of JN . 
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3. Simulations of Quantum Systems 


Quantum systems cannot be simulated with classical computers because 
Hilbert space grows exponentially with the number of particles included in the 
system. Feynman was the first to understand that fact and to propose in 1982 
the construction of quantum computers for the simulation of quantum systems. "? 
Several researchers have since showed that quantum systems can be simulated 


by quantum computers. 


B. PROPOSED REALIZATIONS 


Since 1982 when the field of quantum computation emerged most of the 
research was theoretical and there were no experimental breakthroughs. 
However, in the last few years this trend has changed. Experimental and 
theoretical developments have helped researchers to find and propose some 
possible realizations for quantum computers. A brief overview of two of these 
schemes is presented below. 

1. lon Trap Quantum Computer 

One of the first promising prototypes for quantum computers is the ion 
trap proposed by Ignacio Cirac and Peter Zoller of the University of Innsbruck.”® 
The basis for this scheme is a string of ions stored in a linear ion trap. A 


simplified schematic representation is depicted in Figure 3.1. 
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aubits (ons 


Figure 3.1 lon-trap QC 


In this proposal the internal states of trapped ions represent the qubits. 
The electrostatic repulsive interaction between ions provides the means of 
communication between qubits. Due to Coulomb interaction the linear chain of 
trapped ions has vibrational modes.® The number of normal modes is equal to 
the number of ions. The internal state of each ion can be altered by appropriate 
laser beams, which induce vibrations. These vibrations change the normal 
modes of the system. Thus, we have the transfer of information from one specific 
ion (qubit) to the normal mode of the system. An elaborate treatment of the ion 
trap proposal can be found in the references. 

We will note that there are currently a lot of groups experimenting with 
this realization.® The results from the first experiments are promising.” The 
technology involved is feasible and already being used for other scientific 


purposes.” Finally the most promising aspect of ion traps is that they provide a 
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reliable measurement for the output states utilizing the quantum jump technique. 
This technique can detect the state of an ion with absolute accuracy and almost 
absolute efficiency. 

2. Nuclear Magnetic Resonance (NMR) 

Nuclear magnetic resonance (NMR) is considered today one of the major 
tools in synthetic chemistry.® It is used to study the properties of all states of 
matter except plasma. The field has been developed since 1940's. 

The physical system used for NMR quantum computation® is a single 
molecule with two magnetically distinct nuclei. The qubits are represented by 
nuclear spins, while gate operations are performed by manipulations of the 
spins. The most important method for spin manipulations is the application of 
radio-frequency (RF) pulses. These pulses must have the appropriate power, 
duration and phase. 

The above brief description consists the ideal NMR quantum computer. 
The signal produced from a single nucleus is very small and currently there are 
no available methods to detect it. In order to have a signal strong enough to be 
detected large number of molecules are needed. It is noted that current 
technology cannot detect signals generated by a small number of molecules, 
which is what we need for quantum computing. 

For many years there were two basic problems that made NMR 
computation impossible, the preparation of the initial state and the measurement 


of the output. The preparation problem has to do with the fact that an ensemble 
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of 2-spin molecules in thermal equilibrium contains equally populated spin-up 
and spin-down states. This initial state is totally random and of course is not 
acceptable for a computation where a definite input state must be prepared 
(even as a superposition of other states). The measurement problem arises 
from the fact that for an NMR computer a large number of molecules is used and 
the only output you can measure is the ensemble average magnetization. You 
can't measure the output of a single qubit (molecule). 

In 1997 research groups found solutions for both of these problems and 
since then there is a growing number of groups that search for possible 
implementations.” 

It should be noted that NMR computation has a great advantage due to 
the long coherence times of the spin superposition states. This has as a result a 
reduction in decoherence problems and makes NMR an attractive candidate if 
the other problems will be resolved. 


3. Cavity Quantum Electrodynamics (QED) 

Finally, we note that a very promising approach to develop a universal 
quantum gate exploits the phenomena of cavity QED. In the next section we 
discuss the theoretical background necessary to understand cavity QED. As we 
have noted the operation of a quantum gate is controlled by the Hamiltonian of 


the system. In what follows, therefore, we will concentrate on deriving the 


Hamiltonian for a cavity QED system. 
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IV. CAVITY QUANTUM ELECTRODYNAMICS 


A. GENERAL 


One of the most accurate theories developed in this century is the theory 
of Quantum Electrodynamics (QED), which provides us with the appropriate 
tools to describe and predict the coupling between charged particles and 
electromagnetic fields. The concepts for the interaction of light and matter can 
describe all the physical phenomena except gravity and radioactivity, according 
to Feynman. 

In this chapter we will provide the theoretical background of the 
interaction of atoms with a single mode field. This is necessary for the derivation 
of the Hamiltonian for the system that can be used to implement a 2-bit quantum 


gate. 


B. DERIVATION OF HAMILTONIAN 
1. Quantization of the Electromagnetic Field 


The most elementary approach for the quantization of the electromagnetic 
field is by the means of harmonic oscillator. A field mode is equivalent to a 
harmonic oscillator. This treatment can be found in most quantum optics books 
but we followed the texts by Milonni?' and Loudon.” 

First we begin with a brief review for the harmonic oscillator in quantum 


mechanics.” The Hamiltonian for a harmonic oscillator is given by: 
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2 
P 1 2,2 

H =+—+-—moa 1 
oe ee (1) 


where p and q are operators in a Hilbert space. 


The Heisenberg equations of motion are of the same form as the classical 


Hamilton equations: 


q= Gy '[g,.H =~ (2) 


p=(h)"[p,H]=-mo’q (3) 


Now if we define the non-Hermitian operators a@ and a* 





o- oe (p~ imag) @) 


1 


a= Rear eda (5) 





then q and p can be written as: 


_| R n 
q=i}-*—(a-a ), (6) 
pa @+a") (7) 


and we can write the Hamiltonian as: 





H= sho(aa" +a’ a)j= G +a’ ayho (8) 


The eigenvalues of the operator NV = a*a determine the energy levels of 


the harmonic oscillator, which are given by: 
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E,=(n+3)ho (9) 


The second step is to show that a field mode is equivalent to the 
harmonic oscillator?’. We assume that we work in vacuum where there are no 


sources and the Maxwell’s equations are: 


V-E=0 (10) 

V-B=0 (11) 
1 aB 

VxE=--= 12 

x ae (12) 

vessel (13) 
c Ot 


We define the vector potential A such that B=VxA and the scalar 
potential ¢. Then we can see that equation (11) above is satisfied due to the 


fact that V-(V x A) =0. Also in order to satisfy equation (12) we must have that: 


E=---—_Vo (14) 


V7A-———=0 (15) 
c 


We separate the variables by assuming that A(r,7)=a()A(r). Then 
defining the separation constant k and @=ck and substituting in the above 


equation (15) we get: 
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: sll VA(r) 1%) 2 
a(f)V°A(r) = 2 a(nA(r) > AG) Ta =-k*=> 





> V*A(r)+k7A(r) =0 (16) 


G(t) = -o7 a(t) 


The solutions to the equations (16) give us the following general solution: 
A(r,t) = a(t)A(r) +a" (DA"(r) = a(0)eA(r) + a" (O)e"A*(r) (17) 
Then we can find expressions for the electric and magnetic field vectors in 
the case p=0 (vacuum), using equation (14) and B=VxA, as: 
K(r,1) =-1faaw+a’' an] (17) 
c 
B(r, t) = a(t)V x A(r) + a°(DV x A“(r) (18) 


Then we can find the electromagnetic energy, by its definition, as follows: 


3 2 2 1 lL. 2 3 2 + 7 \2 3A 2 : 2 3 2 
u-ifa HE? +B*)=— (al [ara +4 (0 (ara) +540) farlawf + 


+a(t? [arf x AM +2"? | adrlvx Af +a [av x Aw} 


Now we know that a@(Z) = a(0)e"”, as it was defined above and so: 
a(t) = -iwa(t) a(t’ =-a’ a(t” —s &"(1)? =-w’ a(t)’ 


Also in Appendix A we show that [d°7[V x A@)P =#7 [d*rA(r)’ 
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Finally we can assume that the function A(r) is normalized, so that 
@rA(r)) -1 and the equation for the electromagnetic energy becomes after 
these substitutions: 

a = 1 3 2 2\ _ k? 2 
U=H, =J4 (E27 +B )= leo] (19) 


This last equation is the Hamiltonian for the electromagnetic field. 


Now if we define the quantities: 





_— I es . 

q) == la(r)-a*(2)| (20) 
_— k ba 

PO=T= lair) +a7(] (21) 


the Hamiltonian becomes: 
Hy =>(p* +04") (22) 
Since we have define a(t)= a(0)e" and a*(t)=a"(0)e™ we have that: 
a(t) = -ioa(t), ‘a’ (1) =ioa’ (t) 
From the equations (20) and (21) we can see that g=p and p=-o'q. 


The same result is obtained if we apply the Hamilton equations of motion for the 


Hamiltonian obtained in equation (22): 


gees and peo Sy 
op oq 
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So the defined quantities q and p are canonically conjugate coordinate 
and momentum variables respectively. We can now see the direct mathematical 
equivalence between equation (22), which describes a field mode with frequency 
@ and the equation (1) for the harmonic oscillator of frequency @ and mass m. 


Equation (22) is just a harmonic oscillator with unit mass. 


Now in order to transform the above classical Hamiltonian for the field 
mode to a quantum mechanical one we just follow what we did for the harmonic 
oscillator. Thus we define the lowering and raising operators @ and a™ similar 
to the ones defined in equations (4) and (5), but with unity mass. Then by direct 
comparison with equations (20) and (21) we can find that we have to replace the 


classical variables a(t) and a(t) with respective quantum mechanical operators: 


— 20 - iat} => a(t) = a (23) 
ia*(t) __ + ie ee [2hc? 
ae ire ia’ (t) on >a(hn=a'(t) a (24) 


So we can find that the Hamiltonian for the quantized field mode is: 





#= Sho(aa ee G +atayho (25) 


It is noted that from now on @ and a@* are operators. 
Also the classical vector potential as well the classical electric and 


magnetic fields are now given by new equations and they are now operators in 


the quantum picture: 
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A(r,t) = c ie Jone +a A] (26) 


E(r,1) = i(V2aho JaA(r) - a" (NA) (27) 


2nhe? 





B(r, tf) = Jeov x A(r)+a* (DV x A‘(r)| (28) 


We can also see that the energy eigenvalues are given again by: 
Bz (n+ ho n=0,1,2,... (29) 
We see from this last equation that even for n=0 (no photons) the vacuum 


has energy sho. This is one of the most interesting predictions of QED, the 


existence of zero point electromagnetic field. 

Also we note that in all stationary states the vacuum expectation values 
for E and B fields are zero: 

(E(r, 2)) = (BOr, f)) = 0 (30) 

Up to this point we have talked about a single field mode only. We can 
make the transition now for a multimode field. Even though it is not necessary for 
what we describe later we show here that the proof above can be generalized. 

The conditions we have set so far for the field in free space are that the 
field intensity is independent of position, the spatial part of vector potential A(r) 


satisfies the Helmholtz equation V’A(r)+k’A(r)=0 and the Coulomb gauge 


V-A(r,t)=0. 
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A function that satisfies the above relations is given by: 
A(r) = e"*"é, (31) 
with é, a unit vector and é, -k =0 
If we divide the space into cubes of volume L° and use the boundary 


condition for the vector potential: 








A(x+ZL,y+Z,z+L,t) = A(x, y,2,0) (32) 
we have that : 
2an, 2an 2a, ; 
k, = re k, = ; ~, k= Ft Mas MypM, integers (33) 


In order to satisfy the normalization condition i} @r|A(r)| =1 the vector 


potential A(r) must have the form: 


A, (r)= eee, (34) 


The unit vector é, specifies also the polarization of the field. From the 


fact that é,-k=0 we can see that é, can take two values (with opposite 


directions) é,, ,é,, with: 


A A a a2 a2 _ 
Qn ey =9, en =1 


Each of these values represent the polarizations of the field. So now we 
~ can write the functions for the mode from equation (34) as: 


1 


A, (r) = Hee , (where A=1,2) (35) 
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The equation (26) above for vector potential can now be rewritten as: 


A, (r,t) = ' | 2a Jes (He™ +az,(ne™ i => 


2nhic* 
OV 








>A,(r,0 = Jou (Oe ™ + at, (Oye (36) 


In the above equation @, =kc(as defined before) and a,,,a;, are the 


photon annihilation and creation operators. 
The artificial boundary condition we applied for the field suggests that we 
have an infinite number of modes. Equation (36) describes one mode only (for a 


specific k). So the total potential would be the sum over all possible modes k,, 


due to the linearity of Maxwell’s equations: 





Autd=> 2a Juco raz eh, (37) 


Finally we can write the Hamiltonian as the sum over all possible modes 


and it will have the form of the Hamiltonian for one mode: 
1 r . 
H, = x(5 +A hy jro, (38) 
kh 2 
From our assumptions for the normalization of the modes and the 


definitions for the operators a,,,a;, we can write the following important 


relations: 


fara, Ada = Oe On4! (39) 
v 
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[a,, (1), 02,0] = 534.5, (40) 

In our approach A,, represent an infinite number of discrete mode 
functions that satisfy the artificial imposed boundary conditions. It would be nice 
to have continuous mode functions, but we must point out that this does not add 
anything in the physical meaning. The quantization of the field can be done by a 
different approach with quantum field theory, which is a more elaborated 
treatment and also gives the continuous mode functions. In any case we will not 
need later the quantized field for infinite number of modes but for a single mode 
only. 


2. Interaction of Atoms with Single Mode Field 


Now that we have derived the quantization of the electromagnetic field in 
vacuum we can take the next step which has to do with the derivation of a theory 
that can describe the interactions between atoms and the field. This is the non- 
relativistic QED theory.”" We are still assuming that we work in the Coulomb 
gauge, while we will make some further assumptions in order to reduce the 


terms that will appear in the equations. 


In our derivation we will use again the Hamiltonian and the Hamilton 


equations of motion. It is known that for a charged particle in a field the 


Hamiltonian is: 


H=p-v- mv? +ep—2A-v (41) 
c 


where p we must remember is the canonical momentum described by: 
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p=mv+ <A (42) 
so that the Hamiltonian can take the form: 


ne =(p-£) +ed (43) 
2m c 


Now this last equation determines the time evolution of the particle's 
motion. Our goal is to find a Hamiltonian that will determine also the evolution of 
the electromagnetic field. To do that we assume that the combined Hamiltonian 
should have two terms, one obtained from the electromagnetic energy of the 
field and the other obtained from the kinetic energy of the particle. 

We can see that the kinetic energy of the particle can be written by the 


equation (42) above as: 
1 ee ee | ey 
—mv" =—|p-—A (44) 


We also recall the Hamiltonian found for the field, so we can add these 


terms to express the total Hamiltonian as: 


fe aN 
2m 


(p-£a) +— [a°r(E? +B) (45) 
c 87 

Comparing equations (43) and (45) above we can see that there is a term 
missing, specifically the term eg. In Appendix B we present a treatment, which 
proves that this term is included implicitly in equation (45) and that if we want it 


to appear explicitly in the Hamiltonian we have to express differently the form of 


the electromagnetic energy of the field. 
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As shown in Appendix B the Hamiltonian for the system of particle and the 


field can be expressed as: 


2 
1 e 1 
H=—| p-—A| +e¢+—|d’7\E? +B’ 

mG c ) e an! r( 7 ) ie) 


where E* is the transverse part of the electric field vector. 

We wish now to simplify further the above Hamiltonian using the electric 
dipole approximation.”"” In this approximation we can assume that the vector 
potential A is independent of position and so we can ignore the spatial variations 
of vector potential when we expand the square bracket in equation (46). So we 


can write the Hamiltonian as: 





H=P-£ apr Sia sve@tLfarer +B) (47) 
2m me 2mc? 82 


where we have introduced V = e¢. 
Equation (47) is known as the minimal-coupling Hamiltonian. The 
interaction between the atom and the field is contained in the terms: 


2 
e e ‘ 

Ape A 

mc P 2mc? 





We also note that there are higher order interactions (quadropolar, etc.) 
but their contribution is very small and can be neglected. This is the reason we 


used the dipole approximation for the derivation of equation (47). 
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After all these steps we can now proceed with the quantization of the 
derived classical Hamiltonian. We define the following Hamiltonians for the 


particle and the field respectively: 


H, =P—40(x) (48) 
2m 

Hee ¥ (aie +t ho, (49) 
a 2 


Then we replace the classical variables E, B, p and x by quantum 
mechanical operators. So the total Hamiltonian of equation (47) can now be 


written using equations (26),(48),(49) as: 





H=H,+H, ee (ae Jena eau 





yy me [es 2 Bae te (50) 


"Bme a ka’ 0,0 


It should be noted that since we have assumed the condition of the 
Coulomb gauge we only quantize the transverse electromagnetic field in 
equation (50). The energy of the longitudinal field has been replaced by the 
instantaneous Coulomb interactions among the particles (please see also 
Appendix B). 

In the Hamiltonian of equation (50) we notice that the radiative part is 
quantized and expressed in terms of photon creation and annihilation operators. 


it would be useful to do the same with the terms of Hamiltonian that involve the 
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particle and their interaction with the field. Then we would have the entire 
Hamiltonian in terms of destruction and creation operators for both particle and 
field states. Before we proceed with the quantization we notice that the last term 


; 
A’ in equation (50) can be neglected. This term is already quantized and it 





2mc? 


will not be affected by our second quantization, since it does not explicitly 


involve any atomic operators. Thus, its contribution to the coupling between the 


atom and the field is very small.” 


We can now recall that the completeness theorem® for the basis vectors 
|n) requires that }"|n)(n|=1=0o,, and that H|n)=£,|n), (n|F\m)=F,,,, 
where |”) are the energy eigenstates, F is an operator and o is the density 


operator for particle states. 


Then we can express the Hamiltonian for the particle, applying the 


completeness theorem twice, as: 
at, =| Soo | Dmxnl]=SpXonle Im) = DPN) = DX = Lew 
We can develop a similar expression for the momentum operator p: 


p= | Sool | DinXenl|= DNebnYen| = SrPeanlnl= SrPanl PX |= DPannn 


Thus, we can write the total Hamiltonian as: 





2 
H=%'o,,£, +P (abou +3 po, -=_>y (22% Jl + 55, Pron Fnm (51) 
n kA 2 mc kA nm V : 3 


Q, 
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Now we define the coupling constant for the interaction called g as: 


e | 22 i 
Sanm oe m ho, Pim 7 Cn, (52) 


It is noted that in equation (52) there are no diagonal terms, since they 
vanish due to the fact that the operator pam has odd parity.”* This fact will be 
used below when we will express this Hamiltonian for the two-state atom. 


Then we can rewrite equation (51) as: 
+ 1 + 
H= YO mE, + [atau +s ho, —-hDY YD Bain (c,, +Q, b, (53) 
n ki kA n,m 


This Hamiltonian is now quantized and expressed in terms of operators 
only. Now we can proceed and apply this for a two-state atom. 


The atom has two accessible states |0) and |1). Then equation (53) can 


be simplified. The energy eigenvalues for each state are E) and E; and we have 


that E, -E, =ha,. We also know that o,,+0,, =1, from the definition of o,,, 


above. Then we have: 


1 1 
> £6 mn = E,0,, t Eom = a (o;, +169) +5 Ey(Om +1-9,,)= 


1 1 1 1 
= -Ey\On —Fy)+ oh el ad ao +E,) (54) 


In equation (54) we defined: 


O, =O — Fy (55) 
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Also we have that: 


> Srinm? am = Sxaor%on + Sxav1% 1 = —Sxanr (For — Fro) (56) 


nm 


Now we define: 


O,= (Og, ~ F9) (57) 
0, =F, +O (58) 
Si = Baro (59) 


Then using equations (57) and (59) equation (56) can be written as: 


>, Siam? nm = “Burl y (60) 


It is shown in Appendix C that the two-state operators O,,0,,0, Satisfy 


the Pauli algebra for spin-1/2 particle: 


[o,.0,]= 2io,; Io,,c,]= 2ic,, [o,,0,]= 2io, (61) 
Oo, =0,=0; =1 (62) 


We can also define the operators: 
F =O: O* =O (63) 


Then we have the following properties for these operators: 


oj1)=l0Xil1)=]0) ———@fo)=[1X0l0)=[1) (64) 
oj0)=[0il0)-0 oft) =hYoh)=0 5) 
lo*,o]= Oo, lo,o,]=20 (66) 
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From the above properties, o and o” are called atomic lowering and 


raising operators respectively. 
We are now ready to write the Hamiltonian of equation (53) in terms of 


these operators using equations (54), (56), (57), (60) and (63) as: 


H= sho, +> he, a, Oy, +ih> Sy, (a, +a, (6 -o*) (67) 
kr kA 


In equation (67) we omitted the terms = +£,) and D5hoy,, since 
ki 


there are constants and do not affect the evolution of the system. 
Finally, since we are interested in the interaction of atom with a single 


mode field we can further simplify equation (67) and write it as: 
H= sha,e, +hoa*a+ihg(a+a*\(o-o*) (68) 


This last equation is our final Hamiltonian that describes the interaction of 
the two-state atom with a single mode field. This is our starting point for 


describing the evolution of the proposed cavity QED quantum gate. 
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V. REALIZATION FOR 2-BIT UNIVERSAL GATE 
In the previous Chapter we derived the Hamiltonian for the interaction 
between a two-state atom and a single mode field. In this Chapter we introduce 
some further approximations and we derive the Hamiltonian of the so called 
Jaynes-Cummings model.7" Then we calculate the upper and lower state 
occupation probabilities for the two-state atom. Finally we describe a proposed 


scheme for a universal gate based on the above model. 


A. CALCULATION OF STATE OCCUPATION PROBABILITIES 


As a starting point we use the Hamiltonian of equation (68) of Chapter IV: 
H = Shaye, +hoa*a+iig(a-+a*\o-o") (1) 


One of the most useful approximations in the theory of two-state atoms is 
the rotating wave approximation (RWA).™ This approximation is used here to 
reduce the Hamiltonian of equation (1). The RWA applies to the third term of 
above equation, which is the interaction term between atom and the field. We 
expand this term and find: 

(a+a*\o-o*)=ac-ao*+a’o-a‘o* (2) 

Now we consider the evolution of these operators using the Heisenberg 
picture without the coupling constant g (g=0). 

The time dependence for the photon annihilation and creation operators 


has already been introduced in Chapter IV as: 
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a(t) = a(0je"™ a‘ (t)=a"* (O)e™ (3) 
The atomic lowering and raising operators have a similar time 


dependence: 
o(t) =a(0)e*™ o*(t)=07 (0)e"" (4) 
Then, in the absence of coupling each term of equation (2) above 


becomes, using equations (3) and (4): 


aa* =a(0)o*(0)e (5) 
ato = a* (0)o(0)eo™™ (6) 
aot =a (0)o* (0c (7) 
ao = a(0)a(0)e"* (8) 


We can see that when we are interested for the case when o=a, 


(resonant and near resonant frequencies), the terms @*o* and ao in equations 
(7) and (8) are rapidly varying and their average tends to zero. These terms can 
be neglected. In contrast we see that the rest two terms ao” and a‘o in 
equations (5) and (6) have a different time dependence and their average value 
near resonance is not zero. Thus we keep these two terms in the Hamiltonian. It 
is also noted that we have assumed the coupling between atom and field is very 
weak (g is very small number) otherwise the time dependence of the operators 


would be different and the RWA could not be applied. 
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We can also note the physical meaning?"** of the RWA. The term ac* 
represents the absorption of a photon and the excitation of the atom from lower 
to upper state. The term a*o represents the transition of the atom from upper to 
lower state and the emission of a photon. The term ao represents the 
absorption of a photon and the transition from upper to lower state, while the 
term a*o* describes the emission of a photon and the excitation of the atom. 
There are cases where the RWA can not be applied*’ but for our purposes is 


suitable. So we can now write the Hamiltonian of equation (1) as follows: 
H= Sha, +hoa*a-—ihg(o*a-a‘o) (9) 


This last equation is the Jaynes-Cummings model (first presented by 
Jaynes and Cummings in 1963). Jaynes and Cummings used the dressed-state 
formalism to calculate the occupation probabilities for the atom.2' We know that 
a given problem in quantum mechanics can be treated in different representation 
forms.°”° In the Schrédinger picture the time dependence is carried by the wave 
function, while in Heisenberg picture the time dependence is carried by the 
operators. The form of the Hamiltonian is the same as equation (9) for both 
pictures, as it is showed in textbooks”, except the fact that the operators would 
be time dependent (for Heisenberg). Finally we note that another picture is the 


interaction representation”, where both operators and wave function are time 


dependent. 
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We will work in the Schrédinger picture to calculate the occupation 


probabilities. 


Since we are interested for the interaction of an excited atom with the 


single mode field, we assume that at time t=0 the atom is in the excited state |1) 
and the field has n photons present so it is in the state |”). The form of the wave 


function used is: 


i(nety 


|v) = (4,012) + Buna (| 02 +1))ye 2 (10) 
We use the Schrédinger equation: 


nv) vO) = Hly(s)) (11) 


We will substitute equation (10) into (11) and then we will project on the 


states |1)|”) and |0)|7+1) to get the equations for Ain and Bon. We know that 


* and |B,,.:| will give us the occupation probabilities for the two states. 


the |A 








In 
So first we calculate the left hand of equation (11): 


“iC oe 


in LO) - (4, (of) + Bous()0n-+1)p 


i(ms aot 


7 inn + 5 Jol (I1)|7)+B,.u@On+1)))e  ? (12) 


Then we have to calculate the right hand of equation (11). To do that we 
will use the following relations for the atomic operators who are proved in 


Chapter IV and Appendix C: 











o|1)=|0) — o|0)=0 o*|0)=|1) oo *|1)=0 (13) 


=o) oft) =f) (14) 





Oo, 


Also we know the equations for the field operators: 


aln+1)=n+]|n) a*|n) =Vn+1|n+1) (15) 
Then we can write for the right hand side of equation (11) using (13)-(15): 
~i( meat a 


H\ly@) = (Fe, - rho \( A, (0 \n)e 


~i (ae Sa 


+ | x ho, +(n+ Ihe (B. (A)|0)|n + 1)e + 


; 1 
~i(n+—)ot 
i(n >> 


+(-ingA,,()o*a|1)|n) - ihgB,,,, (No* al 0)|n +1) 


igZls 
+ (nigd,,(t)a*o|1)|n) + iheBy,.,(a*o|0)|n-+ 1) "?™ = 


1 
—-i(n+—)at 
¢ >> 


=> H\y(t)) = (3 ho, + no |, (01) n))e 


cot 


+ - sho, + (n+ ine |e. (0)\n+ ir? + 


: da 
+(ingB,,..(Qvn+1|)|n) + ihgd,,(QVn+1|0)\n+1)e "?" (16) 
Now we substitute equations (12) and (16) in (11) and project onto the 


states |1)|7) and |0)|7+1). Then due to the orthogonality relations we have the 


following equations: 
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1)|n) + Bys(O)0n-+1) sae 


= (n{( (3 ho, +nhe 4,0) I)|n) + [Fhe +(n+ 1h Ba (A)|0)|n +1 aia 7 


in(n|(t [4.0 i) + Boon 0) +1)) —iRGn + a4, (0) 


ing (n{alf(Vn+1B,,.,(0|1)\n))- (Vn +14,, (0}0)|n-+1))p 2 = 


=> inA,, (t) + (n + 5 ha, = -(3 @)+ no }n (0) ~ ihgVn +1B,,.,.0=> 
= 4, () == (0-0) - gn + 1B, (0) (17) 


et (of 4,00 1) + Boga O]O}n-+1)) ~in(n+~Ye0(A, (ft) + Baya (OO) n+ i) ve 
=(n+1(0 (3 ho, + mo \A, (0)|1)|n) + (-3 —ho,+(n+ ine} ones (2) 0)| 72 + ) ra 


~ ihg(n+1[(0lf(V2+18,,.(0)1)\n))—(vn-+14,, (0) 0)|n-+1))p "= 


at 1 1 CGA 
=> ThBy (0) + [n + 5 hen (= G +1)o- 7% Bo (1) + ihgvn+1B,,.() > 


. i ae 
=> Bona 0) = -3@ — By) Bonu (1) + gvn+14,,(0) (18) 
Now since we are interested for the case of resonance where o =a, 


equations (17) and (18) take the final form: 


A,,(t) = ~gV/n + 1Boja (2) (19) 
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By (t)=gvn+1A,@ (20) 
These equations can be solved if we differentiate (19) with respect to time 


and substitute in the value of B,,.,(¢) from (20), and vice versa: 
A,,()=-g°(n+DA4,,( (21) 
Bont) = -8° (2 + I)Bonss O) (22) 
The solutions for equations (21), (22) have the form: 
A, (t) =C, coslgvn +1} +C, sin(gvn+1} (23) 
B,,..(t)=C, coslgvn+1¥+C,sinlgvn+1¥ (24) 
Now if we put the initial conditions for the state of atom-field |1)|7) at time 


t=0 (which means that * =1), we can evaluate the coefficients C;-C4. So we 





A, 





can finally find the occupation probabilities for the two states as: 

\A,,(é)|° =cos*(gVn+1)t [Bona (é)) = sin?(gVn+1)t (25) 

Finally for the case of vacuum inside the cavity (n=0) the probabilities 
become: 


A, (Of =cos’ gt [Boru (| = sin” gt (26) 





The interpretation of equations (26) is very astonishing. They imply that 
an excited atom entering a vacuum cavity (no photons present) can make a 
transition from upper to lower state emitting a photon. This photon can be 


reabsorbed by the atom. There is a sinusoidal exchange of the energy between 
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the atom and the cavity (the so called Rabi oscillations between excited and 
ground state). This model of spontaneous emission is very different from the 
well-known exponential decay in free space. In this case the atom interacts with 
a continuum of modes and not a single mode field. The single mode interaction 
is the phenomenon predicted by Jaynes-Cummings model and is due to the 
vacuum fluctuations”, which are firmly established in the QED theory. The 
emission of the photon is reversible. This reversibility makes this quantum 


system a potential candidate for the construction of a quantum gate. 


B. IMPLEMENTATION OF A 2-BIT UNIVERSAL GATE 
1. Description of the Proposed Apparatus 


In 1995 Tycho Sleator and Harald Weinfurter proposed a 2-bit universal 
gate based on cavity QED techniques.’” We will refer to these authors as SW. 
We are now able with the theoretical background presented earlier to describe 
and discuss this realization of quantum gate. 

The proposed implementation involves microwave cavities, Ramsey 
zones and two-state atoms. The atoms are considered as the carriers of the bits. 


Their two accessible states are the ground and excited state represented by |0) 
and |1) respectively. The cavity is either in vacuum state |0) with zero photons 
present or in state |1) with only one photon present. There are two kinds of 


interaction between cavity and atoms. When the cavity frequencya is equal to 


the frequency @, of atomic transition between excited and ground state, the 
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cavity is “on resonance” and the interaction of cavity with atoms is described by 
the Hamiltonian of equation (9) above. It is noted here that in their paper SW 
present a slightly different Hamiltonian for this interaction. When the cavity 


frequency is not w,, then there are no transitions between excited and ground 


state for the atom inside the cavity. The Ramsey zones can produce a rotation of 
the states of atoms in accordance with the frequency and amplitude of the 
induced field. In the proposed scheme the researchers use Ramsey zones that 
produce a rotation of 1/2 of the spin about the y-axis. The operator for this 
rotation is given by o,, which is defined earlier as o, =i(o,, -o,). 

The interaction of an atom that passes through a Ramsey zone R, the off 
resonant cavity and a Ramsey zone R™ is given according to SW by the 
following transformation to its state: 

U = Ro (e' %2)R (27) 

where 9 is the angle of rotation of the atomic state. 

2. Description of a Control-NOT Operation 

For the implementation of a control-NOT operation using this scheme, two 
atoms are required. The first atom is the control bit while the second one is the 


target bit. A control-NOT operation is characterized by the following truth table: 





Table 5.1 Truth Table for Control-NOT operation 
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When the cavity is on resonance we can transfer quantum states from the 
atoms to the cavity. We showed in previous section that when an excited two- 
state atom passes through the on-resonance cavity and there are no photons 
present (vacuum state) there is a transition of the atom from upper to lower state 
and the emission of a single photon. This photon is stored in the cavity long 
enough to interact with the next atom that passes. If the atom is not excited and 
passes through the on-resonance cavity in the vacuum state, nothing happens. 


So we see that the state of the first atom which is the control bit passes to the 


cavity. The cavity is always in the beginning in vacuum state |0}. If the atom is in 
state |0) the cavity remains in that state. If the atom is in state |1), this state is 
transferred to the cavity as described above and cavity is then in state |1). 


The second atom passes through the Ramsey zones and the off resonant 
cavity. The state of this atom is transformed in accordance with equation (23) 


above. When the cavity is in state |1) the state of the atom undergoes a rotation. 
When the cavity is in state |0) the state of the atom does not change. — 
The control bit leaves the cavity in state [0). When the operation is 


complete the gate is ready to operate again. 


3. Discussion 


The above proposal for a quantum gate is feasible with current 


technology, even though it is far from being miniaturized. The transfer of the 
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state of the control bit to the cavity is predicted by the Jaynes-Cummings model 
as noted earlier. Experiments in 1987 have demonstrated the predictions of the 
theory.” The cavity operated at 21.6 GHz with a temperature of 2.5K. The atoms 
used were rubidium excited by laser radiation to the 63P3,2. state. The transition 
was made inside the cavity to the 61Ds,. state. The flux of the atoms that were 
crossing the cavity was low and the cavity had ample time to “relax.” 

The photon released by an excited atom is stored, as aforementioned, 
long enough inside the cavity, thus allowing interaction with the next passing 
atom. it should be noted that we have not included in our treatment the cavity 
decay, which plays a very important role for the time that the photon will be 
present in the cavity. The cavity decay is one of the decoherence mechanisms 
for this implementation, which will be discussed further in the next Chapter. The 
decoherence problems encountered with this specific realization have not been 
addressed, as far as we know, and they must be studied in order to calculate the 
actual computation time of the gate. 

Another fine point is that the cavity has to operate in very low 
temperatures to avoid the existence of thermal photons. The presence of thermal 
photons in the cavity affects the evolution of the system. There is not a 
sinusoidal exchange of energy between atom and photon and so we now show 
the Rabi oscillations wash out.”4”° 


We know that in thermal equilibrium the probability of a photon to be in 


the nth excited state is given by the Boltzmann factor: 
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ge EnlksT 


ar 
>. evn /kpT 
n=0 


P (28) 


n 


Since the energy takes discrete values according to E, =(n + ho, we 
can write equation (28) as: 
he ur Mh 


é 
P = 29 
her Setter Sg nn 
n=0 n=0 


The denominator is a geometric series with sum: 


2 ~nho, bed —ha, f 1 
e kpT a e "% 0 } = 30 
Se i (0) 


n=0 n=0 


and equation (29) can be written as: 


P= e “7ha(1-6 7) (31) 


So now we can calculate the mean number of photons 77: 


Ld ~ho bed ~nho ha, -ho, bad -ho, 4s 
n => nP, -(1-¢ 7 V5 ne kT -[I-¢ i Vest ard (° ‘| = 
n=0 Ole 7) n=0 


n=0 


~h 
9 Mat 


= (1 - eo vest ee ; CE Oy a ay ae => 
en ee (32) 
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And then the probability of having n photons present when the mean 


number is 7 is given by: 


ni \" nh \ ( nm \'(n+1-A n” 
»,=(=75) (-}-(4) (2 A+] eu ea we 


Now we calculate the expectation value of photons as given by the 








Jaynes-Cummings model, using the wave function introduced in equation (10): 





(n(t)) = (y|a*aly) = (aC \(a|+ Bow (O(n +1\(0 | ‘ 
=> (n(t)) = nA,,| + (2+ DIB 


(4,,(0)1)|72) + Bona (1)|0\\2+1)) )=> 
(34) 


ln On+1 











Using the values calculated for |4,,|° 





and |B,,..{_ from equation (26) : 


=> (n(t)) = neos’(gVn +1)t+(n+1)sin?(gvn +1)t = n+ sin?(gVn+1)t 


Since 7 = )'nP, we can obtain for our model: 
n=0 


(n(t)) = SP n+ sin? (gVn+1)t]= 3 ———_[n + sin?(gVn+1)t]=> 


n=0 = (A+ ae 
ee el oe ee ee 
=> (n(2)) =n + (ae (5) sin (gv n+1)t (35) 


Now we can plot (n(t))—7 to see the effect in the evolution of the Rabi 


oscillations of a thermal filed with a mean number of 7 photons present in the 
cavity. Figure 5.1 below is a plot for #=2, m=10 and m=100. For the time t 


we have introduced the dimensionless quantity 7 = ga, with 7 =1 for all plots. 
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Figure 5.1 plot of (n(t))-7 


We see from the above figure that the Rabi oscillations wash out for a 
thermal field. Also we see that the damping of the oscillations is more rapid 
when the mean number of photons is bigger, i.e. when the temperature is higher. 
That is the reason we need very low temperatures to implement this gate. 

Even though the first results from such an implementation look promising, 
we believe that it is not feasible, at least with current technology, to build a 
quantum network with such gates. The most difficult obstacle is to maintain the 
required coherence between all cavities used for the gates. Another serious 
problem is the fact that all the cavities have to be cooled in very low 
temperature. 

In any case this proposal is a first step towards the construction of 
quantum networks. The basic advantage is that it uses 2-bit gates, which are 


more easily implemented instead of 3-bit. Also the system can be used to carry 


out more experiments. 
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A slightly different approach had been presented”° by a research group of 
University of Innsbruck. Their scheme is also based on cavity QED and the 
Jaynes-Cummings model. They utilize a set of N atoms fixed inside a high Q 
cavity. The interaction of each atom with the single mode field of the cavity 
allows the communication between them. Laser beams are also utilized to 
interact individually with each atom, thus allowing the performance of continuous 
operations between various qubits and the construction of a quantum network. 
The group presented this scheme in a paper published in 1995 and they also 
had some first calculations for the effect of decoherence along with a possible 
solution. 

There are also proposals for 3-bit gates using cavity QED techniques.® 
They are based in the interaction of a single mode field with three-state atoms. 
The Hamiltonian used to describe this interaction with the on resonant cavity is 
the same as the one in equation (2). But the implementation of the gate is 
different from the one described above. It is based on the principle of “dark 
states”. Some researchers believe that this proposal will experience less 


decoherence problems than the 2-bit gate. 
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VI. | DECOHERENCE AND CAVITY QED COMPUTERS 


A. GENERAL 


There are a numerous obstacles that need to be overcome before 
quantum computers can be successfully implemented. We will not deal with 
feasibility questions such as whether the implementation is possible with current 
or future technology, since we believe there are presently many theoretical 
questions that remain unanswered. 

In the previous Chapters we described the elements of quantum 
computers and their possible realization making very naive assumptions. All the 
systems used were considered "ideal". Thus, all the operations were ideal, 
something that will never happen in reality. 

We can categorize the errors encountered in a quantum computer as two 
general types.’* The first type of errors can be called internal errors. They arise 
due to imperfections in the original design of the system and in the initial data. 
The design of the system consists of the quantum system implemented for the 
gates and the measurement techniques for the output. The general belief a few 
years ago that there were no ways to overcome these problems has been 
proved false. Quantum error correction methods along with fault-tolerant design 
of quantum systems have been developed, demonsirating that an actual 


quantum computation could be performed in the future. It is beyond the scope of 
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this thesis to describe and discuss further these fields but we must emphasize 
their importance for the future prospects of quantum computation. 

The second category is the so called external errors. This is the most 
difficult obstacle to overcome and may be the one that will make actual quantum 
computing impossible. This has to do with the fragility of quantum information, 
which is a consequence of the decoherence appearing in the evolution of a 
quantum system. The notion of decoherence and its effect on quantum 
computation is the subject of this chapter. The study of decoherence is very 
difficult. First we will give a general definition of decoherence. Then we will 
present some intuitive considerations in order to understand the phenomenon 
and identify its origin. We will continue presenting some of the possible models 
required to study the phenomenon quantitatively and finally we will discuss the 


effects of decoherence on cavity QED based quantum gates. 


B. DECOHERENCE 
1. General Notion and Meaning of Decoherence 


As was noted earlier (Chapter Il) the reversible unitary evolution of a 
quantum system is necessary for the performance of a successful computation. 
In order to achieve this we must ensure that the evolved system is closed. That 
is the system must not interact with the environment. The interaction with the 
environment has as a result the addition of more degrees of freedom, that are 


not included in the Hamiltonian which describes the closed system. Then the 
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evolution of the system is not the one predicted by the original Hamiltonian and 
any quantum information contained in the system is gradually lost. This process 
is called decoherence. 

The phenomenon of decoherence is-so efficient that is very difficult to 
observe it experimentally. Some superconducting systems were recently 
developed for this purpose and they demonstrated the phenomenon.”” 

As noted by Roland Omnés” the existence of decoherence for 
macroscopic objects is due to their large Hilbert space and their crowded energy 
spectrum. According to perturbation theory the wave function changes under the 
influence of any small perturbation. This change is inversely proportional to the 
difference of the unperturbed energy levels. The closeness of the energy levels 
in a macroscopic object results in very small differences between them, which 
consequently has as an effect an enormous increase of any small perturbation. 
The effect in microscopic objects is not so efficient or so rapidly developed but it 
still exists. Microscopic systems do not have such enormous Hilbert spaces but 
decoherence appears due to their coupling with the environment which is 
considered a system with an enormous Hilbert space. 

There are various models that give quantitative results for certain systems 
and with specific assumptions. However, a general model or theory does not 


exist due to the difficulty of the problem. 
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2. Physical Meaning of Density Matrix Elements 


In Chapter Il we introduced the density operator and its matrix 
representation. We will now present hes meaning of the diagonal and non- 
diagonal matrix elements p,, and p,,, respectively for a statistical mixture of 
states.® 


According to equation (20) of Chapter Il p=> P,P,@). From the 
k 


definition of density operator we can write: 


Pan = >, PrP dan (1) 





Also we have that c* = (u, VW.) and so we can write: 


* = (a, |v ve |u,) (2) 





c, 





Then we can rewrite (1) using (2) and the fact that o, =|y,)(y,| as: 


2 














Prn = > Pt ro (3) 
k 
with Ic, * zero or a positive real number. The p,,, is zero only when ail the 
coefficients |c* * are equal to zero. These coefficients represent the average 








probability of finding the system in the state |w,). The probability is defined as 


average due to the fact that the state is undetermined before the measurement. 


if we perform the same measurement N times with the same initial conditions we 


60 








will find the system Nop,,, times in the state |z,). For this reason p,_ is called the 


u,).° 


For the non-diagonal elements of the matrix we can follow the same 





population of state 


derivation as above. We have that: 


Pam = >, Pr(Pr nm (4) 


Also we can write: 





CnC = (Un [We Xe [Mm (5) 


And we finally have using again the fact that p, =|y,)(y, |: 


Pram = D Pence (6) 
k 





These elements express the interference between the states |w,) and 
|z,), when the state |y,) is a linear superposition of these states. From its 
definition we see that p,,, is the sum of complex numbers. Thus, p,, can be 
zero even if none of the products cic** is zero. If p,, is not zero then we can 


infer that a certain coherence exists between the states |u,) and |u,,). If 9, 


equals zero, however, then this coherence does not exist. That's why the off- 
diagonal elements are called coherences.*® 


3. Reduced Density Matrix 


Now we are going to assume that a quantum system (denoted as system 


A) and the environment (denoted as system E) comprises a new "total" system. 
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As was mentioned in Chapter II the Hilbert space for the total system will be the 
tensor product of the individual Hilbert spaces Ha and He. The density operator p 


of the total system is defined for the total Hilbert space. The trace of this 


operator is defined as: 


Tre= DD (u,(A)(u,(E))o(u,(4)\u,(E)) (7) 


Then we can define a new operator p, that will operate only on the Ha 
space. This operator can be constructed by performing a partial trace on E: 
Pp, =Trzp (8) 
The matrix elements for p, are: 


u,(A)) = ¥. ((u,(A)|(u,(E))Jo(u,-(4))|u,(E))) (9) 


m 





(u,,(A) IP. 
Similarly we can define the operator p, =Tr,p and its matrix elements 


are: 


(u,(E)| | thy (E)) = > ((u,(A) (un ED o(u, (AD) 4 2) (10) 


n 


These two operators are called reduced density operators. Their trace is. 
also equal to 1, like the p operator. From equations (9) and (10) we can also see 
that: 

Trp = Ir,(Tr;P) = Tre(Tr,p) (11) 
The reduced operators allow us to make predictions about the possible 


outcomes of measurements for system A or system E alone. The diagonal and 
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non-diagonal elements of the reduced operator matrix have similar definitions 

and properties as the ones discussed above (equations 3-6) for the operator p. 
This means that the non-diagonal terms of the reduced density operator 

for the system represent the interference between various states of the system. 
In the case where we can identify the interaction Hamiltonian between the 


7729 that the system density 


quantum system and the environment, it is shown 
matrix becomes diagonal (i.e. the non-diagonal terms have vanished) in the 
basis defined by the interaction Hamiltonian. This is the effect of decoherence. 
So if someone wants to find quantitative results for the decoherence of a system 
he must find a way to describe the time evolution of the reduced density matrix 
Ps =Trep. 


However this is not easy. Even though we can describe the evolution of 
the total system (AtE) by equation (18) of Chapter II (nz p(t) =[H), p()]), it is 


very difficult (sometimes almost impossible)° to find a similar equation for the 
reduced density operator of the system A. This is the fact that makes the study 


of decoherence very difficult. 


C. DECOHERENCE IN CAVITY QED BASED GATES 


For the cavity QED based gate that was described in Chapter V there are 
two sources of decoherence. The spontaneous emission by the transition of the 
atom from excited to ground state and the cavity decay (loss of the photon). Both 


these mechanisms are damping mechanisms. The Hamiltonian we derived in the 
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previous Chapters did not include these damping mechanisms. It was an "ideal" 


Hamiltonian. But in an actual gate the coupling with the environment is 
inevitable. Thus, the actual performed computation will be different from what the 
ideal Hamiltonian suggests and may include errors. Thus, we have to find 
quantitative results for the effects of decoherence. This way we will know the 
actual evolution of the gate and we can develop methods to compensate for any 
possible errors. 

As was noted earlier, in order to find quantitative results for decoherence 
effects, we must find an expression for the evolution of the reduced density 
operator of the gate. The most widely used method is by means of the master 
equation.”°”**° There are two different ways to implement the decay of the 
system in master equations. The first one is in the case where the decay of the 
system is not observed (a priori dynamics).”°* The second one is in the case 
where there is implemented a continuous photodection process for the photons 
emitted spontaneously by excited atoms and for photons coming from the cavity 
decay. This is called a posteriori dynamics. 

The basic formalism of quantum theory for such photodetection processes 
can be found in published papers*’** and a few advanced textbooks. The 
implementation of this method in a quantum computer based on cavity QED 
techniques has been presented by a research group of University of Innsbruck.”° 
According to their results it is possible to predict the decoherence effects and 


compensate them, thus making feasible an actual computation. Also they argue 
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that this method can be used to derive quantitative decoherence results for any 
quantum gate based in cavity QED techniques. This is true since the dumping 
mechanisms for all these quantum gates are the same. 

However, if it is implemented in the scheme described in Chapter V, we 
believe that it will give results for a stand-alone quantum gate. A quantum 
computer constructed by such gates will experience more problems since there 
are more possible decoherence mechanisms involved, due to the required 
"wiring" between gates. 

Finally, we note that even if we can't find ways to compensate for the 
decoherence, the study can give us an estimation for the decoherence time 


which in fact is the the time limitation for the performance of the quantum gate. 
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Vil. POSSIBLE MILITARY APPLICATIONS 
It has already been mentioned that quantum computation is a new field 
that has been developed for only a few years. One of the weak points of the field 
_is the difficulty of the notion that quantum mechanics. and computer science can 
be combined. As was noted earlier the two fields have been developed 
independently for over fifty years. There is a significant intellectual gap between 
the theoretical scientists of these fields, which results in the difficulty to 
understand the fundamental aspects of each other. | 

This fact, along with the obstacles encountered so far for the 
implementation of a quantum computer, are the primary reasons for the absence 
of proposed applications. 

The only application that has been proposed until now and could have a 
military interest, as far as we know, is Shor's algorithm. The efficiency of this 
algorithm to perform factorization of large numbers could undermine the security 
of computer systems worldwide. Thus, a government that could build a quantum 
computer to run this algorithm would gain a significant advantage over the rest 
of the world. Quantum computers could play the role of first-class weapons in a 
possible future "information war". 

However, this is not the only military application of quantum computation. 
Someone could think of a variety of applications only by reading the article of 


Richard Feynman’, which was the starting point for quantum computation. 


67 














According to Feynman a classical computer can never simulate a quantum 
system. The only way to accomplish this is by quantum computers. David 
Deutsch had a similar approach to the problem and he also made a very definite 
statement for the capabilities of quantum computer.* They can perfectly simulate 
every finitely realizable physical system. 

Starting from this fact we could think some of the possible systems that 
could be simulated and have a military interest. 

One of the most important applications could be the design and test of 
nuclear weapons. There is already a program called Accelerated Strategic 
Computing Initiative (ASCI) sponsored by the Department of Energy.” The 
primary purpose of this program is to establish a capability for designing and 
testing nuclear weapons solely based in computer simulations. The success of 
this ambitious program is based on the fact that a conventional supercomputer 
could perform all the available calculations. Even though we don't know the 
aspects of the program the theoretical question raised by Feynman and Deutsch 
will remain unanswered until the completion of the project. An alternative 
approach could be the usage of quantum computers for performing the required 
simulations. 

Similarly the quantum computers could perform simulations for the testing 
of future aircraft, missiles and other weapon systems. Even though it is 
premature to discuss the construction of quantum computers, the effects of their 


usage in such applications would have tremendous impact in our society. The 
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cost for the development of new weapons could be reduced dramatically 
(assuming that the cost for the construction of a quantum computer would not be 
very high). A more significant consequence would be the elimination of actual 
tests for the capabilities of these weapons. This would have as a result the 
difficulty to identify the technology level of an enemy by observing the conducted 
tests. Thus, quantum computers could be a new threat in the long run for the 
security of every nation. Fortunately or unfortunately, depending on someone's 
point of view, all these are currently science fiction scenarios and they may stay 
this way forever. 

Another possible application of quantum computers is computer vision. 
Computer vision® is a field of artificial intelligence that describes the 
understanding of the structure and properties of the three dimensional world 
from its two dimensional images. The most important fields of computer vision 
are image meee and pattern recognition. Since images contain an 
enormous amount of information most of the currently developed techniques 
utilize methods for data reduction. This is accomplished by edge detection and 
segmentation algorithms. It is inevitable that in the process some information of 
the original image can be lost. The problem becomes more complicated in the 
case of a passive observer” (an observer that has no control over the image 
acquisition process). Quantum computers could analyze the enormous amount 
of data, in very short times, due to their parallelism capabilities. Most of the 


future unmanned vehicles will rely more and more on computer vision for their 
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navigation and the successful accomplishment of their mission. Thus, we believe 
that a further study is required to determine and analyze the implementation of 


quantum computers and algorithms in computer vision techniques. 


70 








Vill + DISCUSSION-FUTURE PROSPECTS OF QUANTUM COMPUTATION 


In the previous chapters we presented some of the aspects that make 
quantum computers a powerful tool and give them the theoretical ability to solve 
problems that current conventional computers cannot. The theoretical 
advancements in quantum theory and the technological breakthrough in 
quantum devices over the last decade give the chance to actually build some 
fundamental gates that can perform a computation. It is noted that these gates 
are far from being a whole computer system that will operate and solve actual 
problems. 

Also, even if we manage to construct an actual quantum computer, in the 
future it will not be the one that will replace the machines we are using today. 
Furthermore, it won't be a low priced desk-size machine, like the conventional 
ones that almost everyone can buy and use. 

There are many obstacles that need to be overcome before one can 
envision the implementation of quantum computers, the most important being the 
phenomenon of decoherence. For the rest of the problems as are described in 
published papers and articles we have a slightly different opinion. 

We won't argue that difficulties regarding the cost and size of proposed 
machines can force the research community to abandon the effort. The evolution 
of current conventional machines is an example. The continuous development of 


miniaturization techniques transformed the huge, slow processing machines of 
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1950 to today’s mini, ultra fast supercomputers. There is no reason (at least a 
theoretical one) that this could not happen for the proposed schemes of quantum 
computers. 

Even if someone suggests that we could never achieve such 
miniaturization, quantum computers could be built by large companies or 
governments that are interested to solve specific problems, regardless the 
required cost and manpower. 

We believe that quantum computation is an important development and 
its implications in physics and computer science are just beginning to be 
apparent. The field has been emerging for only a few years. Thus, a definite 
assessment for the future prospects cannot be made. In any case we would like 
to point out a few things regarding the current status of the field. 

The absence of proposed applications is a weak point of the field. The 
identification of additional possible applications will attract the investment of 
more resources, which are critical for the continuation of research. 

The proposed quantum algorithms are few and there is an urgent need for 
the development of more algorithms. 

The theoretical research should focus more on the decoherence 
problems encountered in every proposed realization. A possible solution could 


be found by the usage of quantum error correction codes® and by the design of 


fault tolerant quantum gates*. 
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The construction of an actual quantum computer depends on 
technological breakthroughs in quantum devices. However, the miniaturization of 
quantum gates is not a critical problem, at least for the moment, as was already 
argued. 

This thesis was a first approach to describe the fundamental elements of 
quantum computation and their implementation by means of cavity QED 
techniques. It is far from being an elaborate treatment for quantum computers. 
There are as aforementioned more fundamental elements and proposed 
realizations that need to be presented and analyzed. 

Quantum computation is a vast subject involving a lot of research in the 
forefront of theoretical physics. Thus, we believe that a research program should 
be adopted in order to establish a firm result for the feasibility of quantum 
computation. 

The major milestones of such a theoretical program could be the 
following: 

1. Evaluation of all proposed realizations and identification of the most 
promising one. 

2. Further analysis of this realization by identifying the actual computation 
times due to technology limitations and decoherence effects. 

3. Development and implementation of quantum error correction codes to 


achieve the desired computational times. 
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4. Development of appropriate quantum algorithms, which will test the 


computing performance. 

The success of such a program cannot be guaranteed but nevertheless 
we believe its adoption is worthy. 

Even if quantum computation fails, the insight we have already gained for 
quantum systems is very important. We dealt with theoretical aspects that 
maybe we would never encounter or maybe only a few would have. The only fact 
that everybody seems to agree is that we need more years of development to 
reach at a point that we can definitely evaluate the concept of quantum 
computers. And until then we are sure that more profound developments towards 


understanding nature itself will be made. After all that is the ultimate goal of 


physics. 
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APPENDIX A DERIVATION OF RELATION f ar|V x A(r)f =k? @rA(r) 


We want to prove that: 
fa@r[Vx A@P =” f @rA(ry (1) 
We will use two vector identities: 
V-(FxG)=G-VxF-F-VxG (2) 
Vx(VxA)=V(V-A)-V7A (3) 
Then we can write the right hand side of equation (1) using the identity (2) 
with G=VxA and F=A: 
(Vx Ay =(VxA)-(VxA)=V-(AxVxA)+A-Vx(VxA)= 
=V(AxVxA)+AlV-(V-A)-V?2A]=V-(AxVxA)-AV?A (4) 
We will also use the fact that we are working in the Coulomb gauge, in 
which the vector potential A satisfies the following relation: 
V-A=0 (5) 
and Helmholtz’s equation: 
VA+K7A =0 (6) 
Then equation (4) can be written using (5) and (6) as: 


(VxA) =V-(AxVxA)+#7A (7) 
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So now the integral of equation (1) can be written using (7) and the 
divergence theorem as: 

far[v x A@p = [a?rv-(AxVxA)+ [a?rk?A(r) = 

= fdS-fi-Ax(VxA)+ faerk? Ary =” [a?rA@e’ (8) 

The first integral in above equation (8) is a surface integral and vanishes 
due to the periodic boundary conditions imposed on vector potential A. Thus, we 


finally have the required equation: 


farivxA@mf =” [ara 
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APPENDIX B DERIVATION OF EQUATION (46) OF CHAPTER IV 


We start from Helmholtz’s theorem: Any vector field F can be expressed 
into a transverse and longitudinal part with the following properties: 
V-F'(r)=0 VxF’(r)=0 (1) 
Now if we apply this theorem for the E field and we use the fact of the 


Coulomb gauge along with the vector identity VxV¢=0, we can see that the 


electric field vector E = =. Vg, can be decomposed in the two parts as: 
Cc 


E! --.5 and E! =-V¢ (2) 


Then we can write for the E field contribution to the total field energy: 
farre? = (d°r(e? +E’? +2E'E’) (3) 
Now we have that: 
} @ rE: -E’ =0 J 2 rE’? = J ar(V¢) (4) 
Using the vector identity Vi¢dy)=(V¢)w+¢Vy we have that: 


V(GV $) = (VO)(VG)+9V- Vb => (V9) = VGV9)-9V'd (5) 


Now using equation (5) the second equation from (4) can be written as: 


f d°rE'” = (d°r(V9) = [a*rvV(gv¢)- f @rgV'¢ (6) 
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For the first term in equation (6) we can use Gauss theorem to convert the 
volume integral to a surface integral and then let the surface extend to infinity, 


thus making the value of the integral zero: 


[a°rVOV9) = —fdspx =0 (7) 


Now for the second term of equation (6) we have that V’¢ =—4zp and so 


we can write it as: 
An J arog (8) 


The complete Hamiltonian for a system of n charged particles and the 


electromagnetic field can thus be written: 


H= Sate. a.) ge Froese) +s, [art r(E? +B?) (9) 


i=l 
So now since we do not have a continuous charge distribution but 
particles at specific positions we want to find an equivalent expression for the 
second term of equation (9) above. 
We know that the potential @ that one of the charged particles “sees” is 


produced by the rest of the particles (n-1) and can be expressed” as: 





on) =< ro ee (10) 


oF 
Also we know that the product of scalar potential and the charge of a 
particle is the potential energy. So we can write for the potential energy of a 


charge qi: 
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_aF q; 
te Ages Ir, -¥, (11) 


Then the total potential energy of all! particles is: 


ail : 495 (12) 
4m fe -, 





Equation (12) can be written in a different form if we sum over all i, j and 


we divide by a factor of 2 to avoid double counting: 





ms 414; 
dare pa (13) 


ee 
We also know that a discrete set of point charges can be written as a 


charge density using delta functions:”° 


pt) = g.0r-1,) (14) 


i=l 

if we substitute the above expression in the second term of equation (9) 
and use the definition of delta function we get the same expression of equation 
(13). 

Now in the case where all particles are fixed in position except the one 
that we want to know its evolution, then equation (13) represents the potential 
energy V(r) for the moving particle due to the other fixed particles*’’*. This 
potential energy can be written as: 


— V@)=eg) (15) 
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where we have assumed that the particle has charge q=e and the 
potential g(r) is known. Then the Hamiltonian for the particle we are interested 
takes the final form: 


1 e,) 1 
H =—|p-—A| +e¢+—[d*7(E? +B’ 
5 [p : } ep +—| rE +B?) (16) 


This is equation (46) of Chapter IV. 
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APPENDIX C PROPERTIES OF TWO-STATE OPERATORS 


We show that the two-state operators o,,0,,0, satisfy the Pauli algebra 


for spin-1/2 particle. We use the definitions: 


O, = Fo FOr 

Oy =O, — F1) 

O, = 91; — 8% 
Then we have: 
le,.0,]= 0,0, —9,0,, = 1G, + O19 (Og, — Fy) — (Oo, — Fo (Oo, + P19) = 
= i(—Ow +01) —H(Oo — 91) =i0, tio, =2io, (1) 


Similar we can find that: 
b,,0,]=2I0, (2) 
[o., o,J= 2io, (3) 
Also, since we have defined that o,,+0,, =1,we have the following 
properties: 
O2 = (6) + Fp (Oo, +19) =F + Oy, =1 (4) 
Similar we can show that: 
o?=07=1 (5) 


Also we have that: 


o,,|0) = (0, - ©)|0) S 019% oi|0) = 04 0|0) 0190 - 01) = -|0) (6) 
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Similar we can find that: 


o,|1) =(oy ~ Op) 1) = ©19Fo1|1) — Fer1Pr0|1) = 00) — 010 os |1) (7) 


82 








10. 


12. 





LIST OF REFERENCES 


Richard P. Feynman, Int. J. Theor. Physics, 21, 467 (1982). 

Richard P. Feynman, Optical News, February 1982. 

David Deutsch, Proc. R. Soc. London A 400,97 (1985). 

David P. DiVincenzo, Quantum Computation, Science Vol. 270 13 October 
1995. 

P.W. Shor, in Proceedings of the 35th Annual Symposium on the 
Foundations of Computer Science, edited by S. Goldwasser (IEEE 
Computer Society Press, Los Alamitos, CA), p. 124 (1994). 

Claude Cohen-Tannoudgi, Bernard Diu, Franck Laloé, Quantum Mechanics 
John Wiley&Sons, (1997). 

Gerald J. Milburn, The Feynman Processor, Perseus Books (1998). 
Hoi-Kwong Lo, S.Popescu, T. Spiller Introduction to Quantum Computation 
and Information, World Scientific, (1998). 

A. Einstein, B. Podolsky and N. Rosen, Physical Review 47, 777 (1935). 

E. Schrédinger, Proc. Camb. Phil. Soc. 31,555 (1935) 

Charles H. Bennet, Quantum Information and Computation, Physics today, 
October 1995. 

Colin P. Williams, Scott H. Clearwater, Explorations in Quantum Computing 


Springer-Verlag, (1998). 


83 








22. 


23. 


24. 


25. 


26. 


S. Raimes, Many Electron Theory, North Holland Publishing Company, 
(1972). 

D.Deutsch, Proc. R. Soc. London A 425,73 (1989). 

T. Toffoli, Math Systems Theory 14,13 (1981). 

E. Fredkin and T. Toffoli, Int. J. Theor. Physics,21,219 (1982) 

Tycho Sleator and Harald Weinfurter, Physical Review Letters Vol. 74, 
4087 (1995). 

L.K. Grover, Physical Review Letters, 79,325 (1997). 

J.I. Cirac and P.Zoller, Physical Review Letters Vol.74, 4091 (1995). 

C. Monroe et. al. Physical Review Letters Vol.75,4714 (1995). 

Peter W. Milloni The Quantum Vacuum An Introduction to Quantum 
Electrodynamics, Academic Press, (1994). 

Rodney Loudon The Quantum Theory of Light Oxford University Press, 
(1983). 

Stephen Gasiorowicz, Quantum Physics John Wiley&Sons, (1995). 

P. Meystre, M. Sargent III, Elements of Quantum Optics Springer-Verlag, 
(1990). 

Peter W. Milloni and Surendra Singh in Advances in Atomic, Molecular and 
Optical Physics, Vol.28, Academic Press (1990). 


T.Pellizari, S.A. Gardiner, J.I. Cirac and P.Zoller, Physical Review Letters 


Vol.75,3788 (1995). 


84 








2/. 


28. 


29. 


30. 


31. 


32. 


33. 





Roland Omnés, The Interpretation of Quantum Mechanics, Princeton 
Series in Physics, (1994). 

Wojciech H. Zurek Physics Today, October 1991 p.36. 

D. Giulini, E. Joos, C. Kiefer, J. Kupsch, |-O. Stamatescu, H.D. Zeh, 
Decoherence and the Appearance of a Classical World in Quantum 
Theory, Springer-Verlag, (1996). 

A. Barchielli and V.P. Belavkin, Journal of Physics A 24, 1495 (1991). 
Masahito Ueda, Physical Review A, Vol.41, 3875 (1990). 

Masahito Ueda, Nobuyuki Imoto and Tetsuo Ogawa, Physical Review A, 
Vol.41, 3891 (1990). 

Thomas Dean, James Allen, Yiannis Aloimonos, Artificial Intelligence, 


Theory and Practice, Addison-Wesley Publishing Company, (1995). 


85 





86 





INITIAL DISTRIBUTION LIST 


Defence Technical Information COREE ak ech ee 2 
8725 John J. Kingman RD.,STE 0944 
Ft. Belvoir, Virginia 22060-6218 


Dudley KMox Library.........sssccesessesssesesseeenteentssneessterneatenscnnnsnnsnsenstsenessengt tg 
Naval Postgraduate School 

411 Dyer Rd. 

Monterey, California 93943-5101 


Professor William Maler.............cccccsceeseseesseeeresereeerenseeettesseseenenennenennatet sees 
Chairman, Department of Physics, Code PH/Mw 

‘Naval Postgraduate Schoo! 

Monterey, California 93943 


Professor JAMES LUSCOMDE. ...........cscececsecescsetetenteetetstereeeesenenenentnenensessetet ses 
Department of Physics, Code PHI/Lj 

Naval Postgraduate School 

Monterey, California 93943 


Professor Wolfgang BaeT............cscsecessesssreeereerentenetteereeerensentacanenaaansestnst ens 
Department of Computer Science Code CS 

Naval Postgraduate School 

Monterey, California 93943 


CPT Dimitrios GiammakOpoulos...........-..ssessseseseseeesereter eres sserenestesntntnenenensnces 
20 Paparigopoutou Str. 

Athens, 11473 

Greece 


87 





